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A SIMPLE EQUATION FOR CORRELATING TURBULENT 

i 

HEAT TRANSFER TO A GAS 


by R. J. Simoneau and R. C. Hendricks 
Lewis Research Center 

ABSTRACT 


Q5M8 



Many investigators have employed the conventional Dittus-Boelter 
(Nusselt type) equation modified by a wall to bulk temperature ratio to 
correlate turbulent heat-transfer data for gases flowing through heated 
tubes. It has been found that the reported convective data for hydrogen, 
helium, air, and carbon dioxide can be correlated to the same 
as the Nusselt type correlations by use of the equation 



h = k( p v)£* 8 d" 0 * 2 

where K is a constant for any given gas. This equation suggests that 
for the four gases investigated, a convective correlation for a given 
gas need not include the temperature dependence of the thermal and trans- 
port properties. 

Based on the plots of these reported data, the values of K are: 
for hydrogen, 0.0480, for helium, 0.0200, for air, 0.00420, and for car- 
bon dioxide, 0.00385. 

These data covered the conditions of film temperature from 500° to 
3400° R with wall to bulk temperature ratios from 1.1 to 9.9. No appre- 
ciable dissociation was reported in the data. The minimum bulk tempera- 
ture was 150° above the critical temperature. The bulk Reynolds number 
ranged from 5000 to 1,500,000 with heat fluxes up to 3 Btu per second 
per square inch. Tube diameters varied from l/8 to l/2 inch. 





INTRODUCTION 


The question of how temperature-dependent thermal and transport 
properties influence turbulent heat -transfer correlations has been in- 
vestigated extensively in recent years (refs. 1 to 8). Much of this 
work has been directed toward determining the temperature at which to 
evaluate these properties, A large fluid to wall temperature gradient, 
similar to that experienced in rocket cooling channels, will cause a 
correspondingly large change in thermal properties across the boundary 
layer. For this reason investigators have felt that the technique of 
accounting for these properties should have a marked influence on the 
heat-transfer coefficients. 

Bartz (ref. 9), Wolf (ref. 6), and others have suggested that tur- 
bulent heat transfer is dominated by the mass flow rate. In refer- 
ence 9, the relation suggested is 

* * (PV)£ (1 

(Symbols are defined in the appendix. ) A similar approach is presented 
in reference 6 in a plot of h as a function of (pV)^ in figure 37 
of that reference. The nature of this proportionality is explored 
herein. The approach will depart from the conventional approach of 
assuming that this proportionality is best expressed in terms of the 
standard Nusselt type equation. The dominance of mass velocity suggests 
that in correlations based only on measurements of heat flux, mass flow 
rate, bulk temperature, tube diameter, and temperature difference, it 
may not be necessary to account for thermal- and transport -property 
variations across the boundary layer. To determine this, the measured 
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parametric data of some recent heated tube experiments (refs. 1 to 8) 
were reexamined in this report. 

References 1 to 8 represent a series of heated tube investigations 
of turbulent heat transfer to gases at moderate to large values of AT. 
In the regions covered by these experiments the conventional Dittus- 
Boelter (EFusselt type) equation, modified by the wall to bulk tempera- 
ture ratio, has successfully correlated the data. In many cases the 
properties were elevated at film temperature; however, some investi- 
gators reported correlations based on bulk temperatures. The question 
of which temperature to use to evaluate the properties has not been 
fully resolved. This correlation usually takes on one of the two 
following forms: 


NUb - C Ee°- 8 Pr°- 4 (TyT b ) 


Ku f ' 0 Ee £i Pi f' 4 

= C ReO- 8 PrO- 4 (T f /T b ) 


4 * 


( 2 ) 


03) 


For the investigations referenced, a = -0.50 to -0.55 and b = -0.80. 

The effectiveness of a correlation based on the measured parameters alone 
will be determined by comparing it with the reported correlations of the 
same data using equations (2) and (3). 

The data considered for analysis from the references 1 to 8 in- 
cluded all that used by the individual authors for their correlations 
except the data of references 2, 7, and 8. All the runs from refer- 
ence 2 were used, but only one local station was analyzed. The data of 
reference 7 were omitted because they yielded high heat-transfer coef- 
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ficients. In reference 8 the authors reported that the data of the earlier 
report, reference 7, yielded high values of h and would not he included 
in their correlations. The data used for correlation in reference 8 was 
not presented in that report. References 2, 3, and 5 also did not present 
the measured data, therefore, it was obtained by direct communication with 
the authors. Only the gas data of reference 4 was used for correlation. 

The data analyzed covered the overall conditions of film temperature from 
500° to 3400° R with wall to bulk temperature ratios from 1.1 to 9.9. 

No appreciable dissociation was reported in the data. The minimum bulk 
temperature was 150° above the critical temperature. The bulk Reynolds 
number ranged from 5000 to 1,500,000 with heat fluxes up to 3 Btu per 
second per sq. in. Tube diameters varied from l/8 to l/2 inch. A more 
detailed list of the data domains is given in table I. 

CORRELATION TECHNIQUE 

In the heated tube experiments, such as those of references 1 to 6 
heat flux q, mass flow rate a* bulk temperature V and tube diam- 
eter d, are the control parameters. The only other measured parameter, 
wall temperature T w or AD, is dependent on the control parameters. 

If a correlation is to be based on these measured parameters, it then can 
be expressed in the form 

AD = f(q^oi, v*) (4) 

Pressure is not included in equation (4) because there were no apparent 
pressure-level effects in the analysis presented herein. The hydrogen 
data of reference 1, are platted in figure 1 in order to evaluate equa- 
tion (4). Figure 1 is a plot of the relation 



A dependence of h on T^/l^ is evident in figure 1 for T w /T^ > 3 
and also has teen otserved and reported in the references 1 to 6. Data 
from each reference were plotted separately in the same manner as indi- 
cated in figure 1 to determine the dependence of h on the wall to bulk 
temperature ratio. A plot of h against T^/T^ at a given value Of 
(cd/a) 0 * 8 d" 0 ' 2 : for references 1 to 6 is shown in figure 2. The heat- 
transfer coefficient h varies as (Ty/T^) , where a can be selected 
as approximately -0.5. One explanation for the spread and varying 
slope of figure 2 could be effects of the tube length to diameter ratio 
(l/d). It was decided to leave the L/D correction out of this analysis 
for simplicity. 

In figure 3, the abscissa of figure 1 was modified by multiplying 
the right side of equation (5) by V^/T^. Various values of a be- 
tween -0.4 and -0.6 were investigated and appeared to make no signifi- 
cant change in the spread of the data. Figure 3 represents data from 
references 1 to 6 platted with this modification by the wall to bulk 
temperature ratio included. Some of the investigators of references 1 
to 6 reported average heat-transfer coefficients and other reported 
local values. No effort was made to distinguish between them in fig- 
ure 3. There were no apparent differences. Inspection of figure 3 
yields the equation 

h = K(pV)£’ 8 d‘°- 2 


( 6 ) 



where K appears to be a constant unique to each gas. The values of K 
for each of the gases investigated are listed in table II. These values 
were obtained from the lines drawn through the data of figure 3. Equa- 
tion (6) correlates all the data within ±12 percent as shown in table II, 
column 4. 

COMPARISON TO NUSSELT TYPE CORRELATION 
Accuracy 

An indication of the comparative correlation accuracies of equa- 
tions (2) and (6) may be obtained by inspection of table II. Column 1 
represents the deviation of the data as reported by the authors of 
references 1 to 6. These authors used either equation (2) or (3) to 
correlate their data; therefore, the deviations in column 1 represent 
the accuracy with which a Nusselt type equation, such as equation (2), 
will correlate the reported data. Column 4 represents the accuracy 
with which equation (6) will correlate the same data, the data of fig- 
ure 3. Inspection of the hydrogen data, for example, indicates a 
deviation of ±11 percent using equation (6), which compares very favor- 
ably with an average of ±10 percent using the NuflSelt correlation. In- 
spection of figure 3 shows that the spread of the data, column 2 of 
table II, could be decreased and the data would be better fit by using 
a slope other than m = 0.8. This has been observed by others, and 
m » 0.8 is retained for convention and so that equation (6) can be com- 
pared with the reported Nusselt correlations. 
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Relation With Nusselt Equation 

The terms in equation (2) can be regrouped in the following manner: 


h = C 


' 0.6 0.4 
k c 


, \ 0 * 8 -0.2. . .a 

(pVk d (TjT h ) 


The subscript x on the property grouping indicates equations (2) 
and (7) can represent correlations based on either bulk, film, or wall 
temperature properties. The form will remain the same, but the con- 
stant C, or the exponent a, or both will be different for different 
reference conditions (ref. l). When the Nusselt correlation (eq« (2)), 
is written in the manner of equation (7), it is of the same form as 
equation (6). This should be expected since both correlations are based 
on the same measured parameters. Inasmuch as equations (6) and (7) will 
correlate the data to within the same accuracy, a comparison of the two 
equations can be made to yield the relation 


where 


K = Ccp x = const 


k °.6 c 0.4 


Figure 4 indicates that, with the exception of para-hydrogen, for perfect 
and near perfect gases, cp is a monotonically increasing function of tem- 
perature. Since figure 3 indicates that K is a constant for each gas, 
equation (8) requires that C, the conventional Dittus-Boelter coeffi- 
cient, be a function of temperature that decreases in inverse proportion 
to cp. Figure 5 is a plot of C against T for all the hydrogen data 
of reference 1. These C values were computed for bulk, film, and wall 



- 8 - 


reference temperatures for each run using equation (7) with a value of *■ 
a = -0.5. The dashed line represents cp against T for normal hydro- 
gen using the property data equations presented in reference 1. The 
solid line represents C against T based on equation (8), specifi- 
cally. 


= K = 0-0480 

~ " <p(T ) N _ H2 


( 10 ) 


The average deviation of the data about this curve is +14 to -11 percent. 
This average deviation for the C value of each run about the curve of 
equation (10) is also the average deviation about K = 0.0480 for the 
data of reference 1. It is consistent with that reported in column 3, 
of table II for reference 1. 

Inspection of column 3 of table II indicates that the data of each 
investigator shown in figure 3 do not fall symmetrically about the 
lines described by equation (6). This nonsymmetry was inspected for 
film or bulk temperature effect, which might suggest K is also a func- 
tion of temperature. No such effect was apparent in figure 3. The data 
seemed to fall in a random pattern about the mean line without regard 
to temperature. Figure 5, however, suggests a slight bulk -temperature 
effect on K. The temperature effect does not appear to be strong; 
however, K seems to be an increasing function with temperature. From 
the available data, it is not possible to predict a more accurate K 
as a function of temperature than one gets by assuming it constant. 

An inspection of some of the near-critical data of reference 4 also in- 
dicates that K decreases slightly as the bulk temperature approaches 


. * 


the critical value. 



Extrapolation to Other Gases 


Another question of concern is how well can one go from one gas to 
another or extrapolate to unknown gases using equation (6) compared with 
the Nusselt equation. Assume that, knowing K for hydrogen to be 
0.0480, one wishes to find the K values for helium, air, and carbon 
dioxide. A reference temperature at which to perform the calculations 
can be selected from figure 4. This point should be sufficiently re- 
moved from the critical temperature to ensure that all the gases are 
behaving in the same manner, that is, like perfect gases. For these 
gases, a selection of 800° - 1000° F would seem reasonable. From 
equation (8) a Dittus-Boelter type coefficient can be evaluated for 
normal hydrogen at some reference temperature in the following manner: 


C 1000° R 


K n _h 2 
1000° R 


0.0213 


( 11 ) 


Applying equations (8) and (ll) and the data of figure 4 yields the 


estimates of K given in the following table: 


Gas 

Experi- 
mental, 
from 
fig. 3 

Calculated at 800° F 
reference temperature; 
C(T) = 0.0223 
from eq. (ll) 

Calculated at 1000° F 
reference temperature; 
C(T) = 0.0213 
from eq. (ll) 

Propor- 

tionality 

constant, 

K 

Propor- 

tionality 

constant, 

K 

Error 

percent 

Propor- 

tionality 

constant, 

K 

Error, 

percent 

Helium 

0.0200 

0.0196 

-2.0 

0. 0191 

-4.5 

Air 

. 00420 

.00381 

-9.3 

.00383 

-8.8 

Carbon 

. 00385 

. 00348 

-9. 6 

. 00381 

-1.0 

dioxide 









Thus it can he seen that, for the gases investigated, it is possible to * 
est imat e with considerable accuracy the K of a gas based on the ex- 
perimentally known K of another gas. These results tend to support a 
similar extrapolation to gases for which K is not experimentally known. 
At present, however, this extrapolation must be limited to gases whose 
curves of cp against T are of the same general form as those of 
figure 4. 

CONCLUSIONS 

For the gases examined, and probably for any near-perfect gas, 
equation (6) with a specific constant K for each gas will correlate 
turbulent heat-transfer data successfully within currently acceptable 
accuracy limits. Equation (6) is very simple in form and should be 
particularly attractive in design-type parametric studies. 

It is not to be inferred from this study that the properties do 
not influence the actual heat -transfer mechanism. All that is being 
pointed out is that property variations of a given gas do not influence 
a data correlation based on measurements of heat flux, mass flow rate, 
bulk temperature, tube diameter, and temperature difference. More 
sophisticated experiments are needed to assess the influence of property 
variations on heat transfer over a range of temperatures. These experi- 
ments would have to be designed to assess the relation that the basic 
fluid-flow and heat-transfer mechanism have with the thermal and trans- 
port properties. Apparently heat flux, mass flow rate, and temperature 
difference reflect these relations such that they can be correlated 
effectively independent of properties. To determine why this is so would 
seem to be a good experimental and theoretical starting point. 
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, APPENDIX - SYMBOLS 

A area, ft 2 

a,b exponents 

C Dittus-Boelter coefficient, dimensionless 

Cp specific heat, Btu/(lh mass ) (°R) 

d tube diameter, ft 

h heat-transfer coefficient, Btu/(hr) (ft 2 ) (°R)/(ft ) 

K proportionality constant, Btu/(°R) (lb°* °) (hr®* 2 ) (ft°* 2 ) 

mass 

k thermal conductivity, Btu/ (hr ) (ft 2 ) (°R) (ft ) 

L/D ratio of tube length to tube diameter, dimensionless 

m exponent 

Nu Nusselt number, bd/k, dimensionless 

Pr Prandtl number, pc /k, dimensionless 

XT 

q. heat flux, Btu/(hr) (ft 2 ) 

Re Reynolds number, (pV)^d/p x , dimensionless 

p f V-ud 

® e mod modified Reynolds number, — , dimensionless 

df 

T temperature, °R 

Zd? temperature difference, T w - T b , °R 

pV mass velocity, lb mass /(hr ) (ft 2 ) 

p viscosity, lb maSS /(ft ) (hr) 

cp property grouping, k°* 6 c°* 4 /n°* 4 , Btu/(°R) (lb°;® s ) (hr 0 ’ 2 ) (ft 0 * 2 ) 

cn mass flow rate, lb mags /hr 

Subscripts: 
b 


f 


bulk temperature 
film temperature 
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w wall temperature 
x reference temperature 
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TABLE I. - DATA USED FROM INVESTIGATIONS OF REFERENCES 1 TO 8 


Gas 

Film tem- 
perature, 
°R 

Wall to 
bulk tem- 
perature, 
ratio, 

V T b 

Pressure, 
lb/sq. in gage 

Exponents 
a b 

Dittus- 

Boelter 

coeffi- 

cient, 

C 

Refer- 

ences 

Hydrogen, 

helium 

500 - 1350 

1.5 - 9.9 

30 - 1350 

a -0. 55 


a 0. 025 

1 

Hydrogen, 

helium 

800 - 3400 

1.4 - 6.6 

40 - 100 


-0.80 

.021 

2 

Hydrogei* 

helium 

Hydrogen 

500 - 1500 
570 - 900 

1.2 - 3.0 
1.1 - 1.7 

13 200 - 1000 
215 - 640 


■ -.80 

.021 

3 

A 

Air 

580 - 1700 

1.1 - 2.5 

15 - 65 


-.80 

c .023 

5 

Air, 

carbon 

dioxide 

800 - 1300 

1.6 - 2.8 

100 


-.80 

d .037 

6 

Nitroga) 

helium, 

carbon 

dioxide 

1 

350 - 2000 

1.2 - 1.6 

up to 1000 

-.50 

1 


e . 036 

00 

•s 

r- 


a More than one value of a and C was used in ref. 1. 
^Information received in a private communication. 

°Wall temperature correlation was used instead of film correlation. 
d Wall temperature correlation and L/D correction were used. 

^ube length to diameter ratio correction was used. 
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TABLE II. - COMPARISON OF EQUATION (6) WITH NUSSEIT CORRELATION 


Gas 

Proport ion- 
' ality 
constant, 

K 

Average Deviation, 
percent 

Refer- 

ence 

. 1- 
(a) 

■EMI 

R5I 

mmm 


Hydrogen 

0.0480 

43 

±12 

+13 to -11 

SHS 

• 1 



±10 

±12 

+18 to -6 

m 

2 



f ±5 

::±7 

+5 to -9 

-ii 

3 



— 

±15 

+9 to -20 


4 

Helium 

0.0200 

±8 

±9 

+12 to -5 

, ip -f p. 

1 



±10 

±13 

+17 to -8 

-r-Lo o V 

2 



f ±8 

±10 

+7 to -12 

-10 

3 



±12 

(e) 

(e) 

(e) 

7,8 • 

Air 

0.00420 


±14 

+16 to -12 

+11 to 

5 



f ±7 

±8 

+8 to -8 

-10 

6 

Carbon 

0.00385 



+8 to -9 

+8 to 

6 

dioxide 





—9 


• 

(e) 



(e) 

(e) 

7,8 

Nitrogen 

(e) 

±12 

(e) 

(e) 

(e) 

7,8 


a Deviation as reported in the references based on eq.. (2) 
or (3). 

^Spread of the data of each reference plotted as in fig. 3. 

c Deviation for the data of each reference from the line 
described by eq.. (6). 

^Deviation for all the data of a given gas from the line 
described by eq.. (6). 

insufficient data published to determine K or to make 
deviation analysis. These data appear to support the 
rest in about the same manner; however, they seem to in- 
dicate higher K values than other reported data. (See 
INTRODUCTION) . 

•‘•Estimated from published curves. 
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EFFECT OF MASS VELOCITY 

AND TUBE GEOMETRY ON HEAT- TRANSFER COEFFICIENT 



FIGURE 1. 



EFFECT OF WALL-TO-BULK TEMPERATURE RATIO 
ON HEAT- TRANSFER COEFFICIENT 



FIGURE 2. 




CS-31126 





EFFECT OF TEMPERATURE ON PROPERTY GROUP 
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FIGURE 4. 


EFFECT OF TEMPERATURE ON DITTUS-BOELTER COEFFICIENT 



FIGURE 5. 



